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STATEMENT  OF  THE  PROBLEM  STUDIED 

The  investigations  performed  during  the  granted  period  were  aimed  at 
understanding  of  the  processes  of  inhibition  of  combustion  of  methane  and  propane 
flames  by  organophosphorus  compounds  (OPC),  which  are  known  to  be  perspective 
ecologically  safety  fire  suppressants.  Molecular  beam  mass  spectrometry  (MBMS)  and 
computer  modeling  were  shown  to  the  most  effective  technique  to  study  the  structure 
of  CH4/02/Ar  and  C3H8/02/Ar  flames  doped  with  trimethylphosphate  (TMP)  and 
stabilized  on  a  flat  burner  at  1  atm.  The  foreground  task  of  experiment  was  to  measure 
profiles  of  temperature  and  concentration  of  stable  and  labile  species:  free  atoms  and 
radicals  and  phosphorus-containing  products  of  TMP  destruction  in  flames  PO,  P02, 
HOPO,  H0P02  and  OP(OH)3,  which  are  responsible  for  inhibition  process.  OPC 
loading  dependence  of  burning  velocity  of  atmospheric  CH4/air  and  C3H8/air  flames 
was  planned  to  be  measured  using  Mache-Hebra  nozzle  burner  and  total  area  method 
of  a  flame  image.  Comparing  the  experimental  data  with  modeling  results  and 
applying  sensitivity  analysis  of  the  kinetic  model  for  TMP  destruction  elaborated 
earlier  was  planned  to  be  refined.  Finally  the  refined  kinetic  model  must  predict  both 
burning  velocity  and  chemical  structure  of  an  atmospheric  flame.  Fire  suppression 
effectiveness  of  OPC  for  premixed  methane  flame  stabilized  on  opposed-jet  burner 
was  planned  to  be  studied. 
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'Scientific  Progress  and  Accomplishments 
INTRODUCTION 

Intensive  study  of  destruction  chemistry  of  organophosphorus  compounds  (OPC) 
began  about  10  years  ago  and  was  associated  mainly  with  a  problem  of  chemical  warfare 
agents  (CWA)  disposal  as  incineration  was  recognized  to  be  the  basic  technology  for  CWA 
disposal  in  USA.  The  first  results  were  obtained  in  USA  and  Russia  practically 
simultaneously  and  were  devoted  to  the  study  of  structure  of  Efi/CVAr  flames  doped  with 
dimethyl  methylphosphonate  (DMMP)  stabilized  on  a  flat  burner  at  a  pressure  of  about  50 
Torr  using  molecular  beam  mass  spectrometry  (MBMS)  with  soft  ionization  by  electron 
impact  [1,2]  and  vacuum  ultraviolet  [3].  Kinetic  model  for  DMMP  destruction  in  a  flame  was 
proposed  by  Wemer  and  Cool  [3]  and  was  based  on  qualitative  results  on  structure  of 
lU/CU/Ar  flame  doped  with  DMMP  [3]  and  thermochemical  and  kinetic  data  of  C.  Melius 
calculations  [4]  for  the  stages  involving  DMMP  and  phosphorus-containing  intermediate 
products  of  it  destruction.  To  describe  the  interaction  of  final  phosphorus-containing  products 
with  H  and  OH  radicals  Twarowski  mechanism  [5,6,7]  was  applied.  Comprehensive  study  of 
H2/02/Ar  [8-11]  flames  doped  with  OPC  (trimethyl  phosphate  (TMP)  and  DMMP)  and 
stabilized  on  flat  burner  at  low  pressure  was  performed  by  the  authors  and  resulted  in 
development  of  kinetic  mechanism  [10,12].  The  model  [10,12]  adequately  described 
experimental  data  on  concentration  profiles  of  stable  and  labile  flame  species  including 
radicals  H,  O,  OH,  PO,  PO2  and  other  phosphorus-containing  species  (PCS)  HOPO,  HOPO2, 
OP(OH)3  at  low  pressure. 

Recently  after  signing  Montreal  Protocol  forbidding  production  of  Halon  1301  and 
search  of  candidates  for  its  replacement  OPC  attracted  attention  as  fire  suppressants  and  the 
need  of  study  of  mechanism  of  action  of  OPC  on  combustion  has  emerged.  In  connection  with 
this  problem  the  mechanism  of  action  of  OPC  additives  on  low  pressure  and  atmospheric 
hydrocarbon  flames  was  studied.  Structure  of  methane  premixed  laminar  flat  flames  was 
studied  using  MBMS  and  computer  modeling  [12-15].  Methane  opposed-jet  flames  were 
studied  using  laser  induced  fluorescence  and  computer  modeling  [16].  Effectiveness  of 
number  of  OPC  as  fire  suppressants  was  detennined  using  cup-burner  technique  [17].  The 
investigations  [14,15]  revealed  the  drawbacks  of  early  proposed  mechanism  of  OPC  action 
and  all  efforts  were  exerted  for  its  refinement.  The  investigations  of  Glaude  et  al  [18,19]  were 
devoted  to  the  development  of  mechanisms  of  OPC  destruction  in  flow  reactor  and  in  flames. 
Comparison  of  calculated  data  obtained  with  experimental  results  on  structure  of  low  pressure 
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HyC^/Ar  flame  doped  with  DMMP  and  TMP  demonstrated  satisfactory  agreement  not  only 
for  stable  flame  species  but  also  for  labile  PCS  such  as  PO,  PO2,  HOPO  and  HOPO2  [18]. 

Babushok  and  Tsang  [20,21]  proposed  another  kinetic  mechanism  for  DMMP 
destruction  in  CHVair  atmospheric  flame.  This  mechanism  was  applied  for  calculation  of 
burning  velocities  of  stoichiometric  CHVair  flame  at  atmospheric  pressure. 

In  spite  of  significant  progress  in  study  of  combustion  chemistry  of  OPC  there  is  a 
deficit  of  quantitative  experimental  data  on  flame  structure,  which  are  needed  for  further 
refinement  of  a  mechanism  of  action  of  OPC  on  combustion.  Besides,  the  above-mentioned 
models  were  applied  to  simulate  only  a  several  experimental  data  on  flame  structure  and 
burning  velocity.  The  goal  of  present  research  was  to  improve  our  understanding  of  OPC 
combustion  chemistry  and  the  inhibition  mechanism  through  experimental  and  modeling 
study  of  structure  and  velocity  of  propagation  of  atmospheric  CjHg-TVAr  flames  doped  with 
TMP. 

EXPERIMENTAL 
Experimental  Conditions 

To  investigate  OPC  effect  on  the  structure  of  atmospheric  flames  four  following  flames  have 
been  chosen: 


Table  1 


COMPOSITION  (in  mole  fractions) 

Equivalence 

ch4 

02 

Ar 

ratio  (cp) 

0.06 

0.15 

0.79 

0.9 

0.075 

0.125 

0.80 

1.2 

c3h8 

02 

Ar 

0.025 

0.136 

0.839 

0.9 

0.029 

0.121 

0.85 

1.2 

The  flames  were  stabilized  on  a  flat  Botha-Spolding  type  burner  16  mm  in  diameter  at 
atmospheric  pressure  and  temperature  370-400  K.  Combustible  mixture  was  prepared  with  a 
help  of  mass  flow  controllers  (MKS  Instruments  Inc.,  model  1299S).  Total  volumetric  flow 
rate  was  1 .5  slpm.  OPC  was  added  to  a  combustible  mixture  flow  using  saturator  with  liquid 
TMP  in  controlled  temperature  bath. 

Methodic  problems  Development  of  MBMS  Technique  for  Atmospheric  Flames. 

During  experimental  study  it  has  been  found  out  that  the  time  of  exploitation  of  quartz 
probes  is  much  shorter  in  atmospheric  flames  than  in  low  pressure  flames.  It  occurs  because 


9 


of  2  reasons.  The  first  one  is  explained  by  higher  temperature  of  a  probe  in  atmospheric  flame 
due  to  more  effective  thermal  contact  between  a  flame  and  a  probe  at  atmospheric  pressure. 
When  quartz  is  heated  above  1200°C  and  kept  at  the  temperature  for  the  prolonged  time  it 
starts  crystallizing.  This  process  is  accelerated  by  contamination  of  the  surface  and  is 
followed  a  quartz  probe  destruction.  Besides,  while  sampling  the  phosphorus  compounds 
form  a  solid  film  of  pyrophosphoric  acids  (especially  near  the  tip  of  the  probe)  which  alloys 
with  heated  quartz  to  fonn  easily  melted  phosphate  glass. 

The  systematic  research  of  chemical  structure  of  gaseous  flames  at  1  atm  using 
MBMS  was  not  performed  in  the  past.  It  may  be  explained  mainly  by  methodical  difficulties 
in  taking  into  account  the  appreciable  perturbations  of  the  flames  caused  by  a  sonic  probe 
When  studying  the  structure  of  atmospheric  flat  flames  using  molecular  beam  mass 
spectrometry  (MBMS)  it  has  been  discovered  that  measured  concentration  profiles  of  main 
flame  species  CH4,  O2  and  fFO  are  displaced  from  the  temperature  profile.  These  profiles  are. 
shifted  in  accordance  with  the  shift  of  lines  of  equal  concentrations  [22]  in  a  flame  caused  by 
a  probe  It  is  explained  by  appreciable  heat  losses  from  the  flame  not  only  into  the  burner  but 
into  a  probe  and  surrounding  space.  The  flame  heat  losses  into  the  probe  depend  on  the 
probe's  position  relatively  the  burner  and  can  not  be  estimated. 

To  measure  temperature  profile  of  the  perturbed  flame  correctly  one  should  use  a 
thermocouple  placed  near  the  probe's  tip.  This  technique  was  proposed  by  two  groups  of 
authors  [23,24]  and  was  applied  for  flat  flames  stabilized  at  low  pressure.  Earlier  [25]  we  have 
demonstrated  that  for  the  probes  having  the  wall  thickness  0.05  mm  the  temperature  profiles 
in  low  pressure  flames  measured  by  a  thermocouple  near  the  probe's  tip  and  far  off  probe 
practically  do  not  differ  from  one  another.  The  authors  [23,24]  recommend  to  a  place  a 
thermocouple  at  the  distance  of  2-3  diameter  of  the  thermocouple's  junction  from  the  probe 
that  is  0.2-0. 3  mm  in  the  case  of  the  thermocouple  0.1  mm  in  diameter.  At  atmospheric 
pressure  when  the  effect  of  a  probe  on  the  flame  is  rather  stronger  the  selection  of  optimal 
distance  between  a  probe  and  a  thermocouple  was  performed  by  comparison  of  concentration 
profiles  of  CH4,  O2  and  H2O  with  profiles  of  temperature  measured  by  a  thermocouple  which 
junction  was  placed  at  varied  distances  8  from  the  probe's  tip.  Figures  1,  2  and  3  demonstrate 
a  perturbation  of  the  flat  atmospheric  methane  lean  and  rich  flames  with  additive  and  without 
its  by  a  sonic  probe.  Comparing  the  temperature  profiles  (Figs.  1-2)  obtained  with  a  help  of 
thermocouple  placed  at  various  distances  from  the  tip  of  the  probe  one  can  see  that  the  probe 
decreases  final  flame  temperature  for  about  200  K  and  noticeably  increases  the  width  of  a 
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combustion  zone.  It  is  noteworthy  that  the  concentration  profiles  should  be  moved  towards 
the  burner  on  the  value  of  sampling  shift  [22].  The  goal  of  the  comparison  is  to  achieve  the 
similarity  of  profiles  of  concentrations  of  stable  flame  species  and  temperature.  The 
experiments  demonstrated  optimal  distance  to  be  8=0.3±0.03  mm  for  above  flames  and  probe 
used.  It  is  evidently  that  this  recommendation  can  be  applied  for  certain  types  of  probes  and 
flame. 

Temperature  Measurements 

The  temperature  of  the  flames  was  measured  by  a  Pt  -  Pt+10  %  Rh  thermocouples 
welded  from  wire  of  0.02  mm  in  diameter.  The  thermocouple  covered  by  a  layer  of  SiC>2  was 
0.03  mm  in  diameter.  The  thermocouple  had  a  Il-shape  with  shoulders  about  2  mm  long.  At 
such  ratio  wire  diameter/shoulder  length  the  heat  losses  in  the  cold  ends  of  a  thermocouple  are 
negligibly  small. 

To  verify  the  accuracy  of  measurement  of  temperature  profile  and  estimation  of  the 
sampling  shift  the  method  of  pneumatic  probe  [26]  was  applied  to  measure  the  temperature 
profiles  in  the  flames.  The  temperature  profile  obtained  using  this  technique  satisfactorily 
coincided  with  that  measured  by  a  thermocouple  placed  at  the  distance  of  0.3  mm  from  the 
probe’s  tip  for  above  mentioned  flames  and  probe.  But  for  another  probe  the  temperature 
profile  obtained  using  the  method  of  pneumatic  probe  satisfactorily  coincided  with  that 
measured  by  a  thennocouple  placed  at  the  distance  of  0.6-0. 7  mm  from  the  probe's  tip  (Fig  4). 
Actually  a  pneumatic  probe  detennines  the  averaged  temperature  in  the  sampling  area. 
Therefore,  the  congruence  of  temperature  profiles  measured  using  both  techniques  confirms  a 
correctness  of  the  distance  between  a  thermocouple  and  the  probe,  which  was  chosen  for 
measuring  flame  temperature.  Besides  it  is  necessary  to  note  that  modeling  data  justified  the 
similarity  of  profiles  of  temperature  and  species  concentrations  in  flame.  Accordingly  to  the 
data  of  modeling  of  structure  of  C^Ffi/CF/Ar  flame  stabilized  on  the  flat  burner,  concentration 
profile  of  FFO  is  similar  to  temperature  profile  (Fig.  5).  Thus,  regularity  of  the  choice  of  the 
"proper"  temperature  profile  was  confirmed  also  by  the  similarity  of  the  above  temperature 
profile  and  measured  concentration  profile  of  O2. 

Experimental  Setup 

The  quadrupole  mass  spectrometer  with  a  molecular  beam  sampling  system  has 
been  partially  described  previously  [2,8]  For  the  current  experiments,  the  geometry 
of  the  sampling  system  was  as  follows:  probe  -  skimmer  distance  -  20  mm,  probe  - 
collimator  distance  -  320  mm,  probe  -  ion  source  distance  -  390  mm.  Stage  1  was 
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pumped  by  an  oil  diffusion  pump  (2000  I/s),  which  provided  a  working  pressure  of 
5x-10'4  Torr.  Stages  2  and  3  were  pumped  with  a  turbomolecular  pump  (500  I/s), 
which  maintained  pressures  of  about  10  5  Torr  (stage  2)  and  2xlCf8  Torr  (stage  3) 
during  the  experiments.  The  M  S  7302  quadrupole  mass  spectrometer,  produced  by 
the  Experimental  Plant  of  the  Russian  Academy  of  Sciences,  was  equipped  with  a 
modernized  ion  source  with  a  small  spread  of  electron  energy,  corresponding  only 
to  the  thermal  spread  [27,28].  In  practice,  ionization  potentials  were  measured 
e.g.,  of  argon  and  oxygen  with  an  accuracy  of  ±  0.2  eV.  This  ion  source  allowed 
operation  at  low  ionization  energies  (IE)  close  to  the  ionization  potentials  of  the 
atoms,  radicals  and  molecules  investigated.  The  molecular  beam  was  chopped 
(frequency  33-49  Hz)  and  a  photodiode  signal  was  used  as  a  base  signal  in  the 
measurements  of  the  modulated  component  of  the  molecular  beam.  The 
measurements  were  made  under  ion  counting  conditions  by  CAM  AC  apparatus 
and  a  computer.  The  power  supply  for  the  chopper  motor  was  computer- 
controlled  and  the  rotation  frequency  was  defined  with  a  timer.  The  photodiode 
signal  passed  through  a  phasor,  which  initiated  the  time  interval  for  the  gated  pulse 
counter  to  operate.  The  pulsed  signal  from  the  electron  multiplier  was  amplified 
and  then  passed  to  the  pulse  counter.  The  base  frequency  of  1  M  Hz  for  the  timers 
was  produced  by  a  generator.  The  ionization  energy  of  electrons  was  controlled  by 
a  digital -to- analog  converter  with  a  stability  of  ±  0.005  V.  To  determine  the 
reliability  of  the  measured  concentrations  and  their  reproducibility,  each  one  was 
deduced  from  4  to  6  measurements.  The  measurements  usually  did  not  last  more 
than  12  s.  A  mean-square  error  was  estimated;  only  those  peak  intensities  with 
mean- square  errors  of  less  than  35%  were  used. 

Burning  Velocity  Measurement 

Study  of  OPC  additives  on  burning  velocity  of  flames  was  performed  using  Mache- 
Hebra  nozzle  burner  described  elsewhere  [29]  and  the  total  area  method  from  a  flame  image. 
Burning  velocity  was  estimated  as  u=W/S,  where  W  -  volumetric  flow  rate  of  the  combustible 
mixture,  S  -  area  of  the  flame  cone.  Accuracy  of  measurement  of  burning  velocity  was  about 
4%.  Experimental  technique  was  described  earlier  [15]. The  size  of  the  flame  cone  was 
measured  by  luminous  zone  but  not  with  a  help  of  schlieren  method.  This  technique 
introduced  an  error  in  absolute  values  of  burning  velocity.  That  is  why  here  we  present  the 
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ratio  (f)  of  burning  velocity  of  inhibited  flame  (ug)  to  burning  velocity  of  undoped  flame  (uo) 
as  a  value  characterizing  the  effectiveness  of  an  inhibitor.  To  ascertain  in  reliability  of  the 
data  obtained  magnitude  of  f  was  measured  in  different  way.  For  that  the  height  of  the  cone  of 
doped  flame  was  decreased  and  became  equal  the  cone  of  a  undoped  flame  owing  to  decrease 
of  volumetric  flow  rate  of  the  combustible  mixture.  In  this  case  effectiveness  of  an  inhibitor 
was  estimated  as  f=Wg/W0,  where  Wg  and  Wo  are  the  volumetric  flow  rate  of  combustible 
mixture  for  inhibited  and  for  uninhibited  flames  correspondingly.  Results  obtained  by  both 
methods  demonstrated  satisfactory  coincidence. 

MODELING 

PREMIX  and  CHEMKIN-II  codes  [30,31]  were  used  for  simulation  of 
concentration  profiles  and  burning  velocity.  Simulation  of  flat  flame  structure  was  performed 
using  specified  experimentally  measured  profiles  of  temperature.  A  number  of  kinetic  models 
for  oxidation  of  methane  are  available  (see  Table  2).  The  simplest  one  involves  18  species  in 
59  stages.  This  mechanism  involves  only  Ci  hydrocarbons  and  is  incorporated  in  PREMIX 
and  CHEMKIN  II  software  package  [30,31].  Besides,  there  are  2  models  [32,33]  for  methane 
oxidation,  which  involves  C1-C3  species. 

Table  2.  Brief  Characteristics  of  Mechanisms  for  Hydrocarbon  Oxidation. 

Mechanism  Number  of  C-containing  species  Number  of 


stages  species  N-containing  species 


CHEMKIN 

58 

18 

Ci 

0 

[30,31] 

GRI  3.0  [32] 

325 

53 

C1-C3 

8 

Leeds [33] 

175 

37 

C1-C4 

0 

Konnov  [34] 

1200 

127 

C1-C4,  c6 

35 

For  study  of  C3Hg/02/Ar  flames  at  atmospheric  pressure  the  mechanism  based  on  Konnov's 
detailed  reaction  model  [34]  for  combustion  of  small  hydrocarbons  was  used.  This 
mechanism  involves  127  species  containing  H,  O,  C,  N  elements  and  1200  elementary 
reactions.  Mechanisms  for  destruction  of  TMP  and  DMMP  in  ^/CF/Ar  and  CFLj/CF/Ar 
flames  stabilized  on  a  flat  burner  at  low  pressure  of  47-76  Torr  was  elaborated  earlier 
[10,12,14].  The  model  comprises  14  stages  for  destruction  of  TMP  and  organophosphorus 
intermediates,  5  stages  for  destruction  of  orthophosphoric  acid,  and  15  stages  of  so-called 
reduced  Twarowski  mechanism  [5,6,7],  which  were  modified  by  us  [10,12].  Later  these 
models  were  applied  to  simulate  the  structure  and  burning  velocity  of  atmospheric  flames.  To 
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obtain  a  satisfactory  agreement  between  experimental  and  modeling  results  a  number  of  rate 
constants  were  modified. 

One  of  the  kinetic  models  for  OPC  destruction  in  a  flame  describing  inhibition  effect 
of  OPC  on  combustion,  which  was  used  to  simulate  the  flame  structure,  was  the  mechanism 
proposed  by  Glaude  et  al  [18,19].  A  kinetic  model  of  DMMP  and  TMP  destruction  in  a  flame 
includes  202  reactions  involving  41  phosphorus-containing  species  (PCS).  The  authors 
evaluated  rate  constants  of  many  reactions.  Rate  constants  of  the  reactions  from  Twarowski 
mechanism  have  been  evaluated  too.  Another  mechanism  proposed  by  Babushok  and  Tsang 
[20,21  comprises  24  stages  with  participation  of  OPC  from  the  model  of  Werner  and  Cool  [3], 
79  reactions  from  Twarowski  mechanism  (including  reactions  of  phosphine  oxidation  with 
participation  of  PH3,  PTC,  PH)  and  number  of  stages  of  interaction  of  CH3  and  CH3O  with 
phosphorus  oxides.  It  should  be  noticed  that  our  kinetic  model  for  OPC  destruction  as  well  as 
the  model  proposed  by  Babushok  and  Tsang  are  based  mainly  on  thermochemistry  estimated 
by  C.  Melius  [4].  The  model  elaborated  by  Glaude  et  al.  has  thermochemical  data  estimated 
by  the  authors. 

Structure  of  CH JO  7/Ar  Flames 

Figure  6  demonstrated  experimental  temperature  profiles  in  the  lean  CH4/02/Ar  flame 
doped  with  350  ppm  and  without  additive.  The  measurements  were  perfonned  using  a 
thermocouple  placed  away  from  the  probe.  As  can  be  seen,  in  the  presence  of  TMP  the  width 
of  the  combustion  zone  increases  (i.e.  the  additive  inhibits  combustion)  and  the  final 
temperature  in  the  post-flame  zone  is  slightly  higher.  Visual  observations  demonstrate  that  the 
flame  turns  orange  and  a  distance  between  the  burner  and  luminous  front  increases. 

Thennocouple  measurements  showed  that  the  profiles  obtained  by  a  thennocouple 
placed  at  the  distance  0. 1-0.3  mm  from  the  tip  of  a  microprobe  are  almost  identical  to  those  of 
unperturbed  flame.  It  indicates  that  the  microprobe  introduces  minor  thermal  perturbations 
into  the  flame. 

The  temperature  of  the  burner  surface  was  additionally  measured  by  a  thennocouple 
embedded  into  the  center  of  the  burner's  filter.  The  temperature  of  the  burner  surface  was 
demonstrated  to  be  the  same  for  both  undoped  flames  -  nearly  400K.  For  the  doped  flames 
where  the  preflame  zone  is  increased  in  comparison  with  uninhibited  flames  the  temperature 
is  370K.  The  influence  of  a  probe  on  the  temperature  of  the  surface  is  insignificant  (within  2- 
3K). 
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Figure  7  shows  the  effect  of  350  ppm  of  TMP  on  concentration  profiles  of  CH4,  O2 
and  H2O  in  the  lean  CFfi/C^/Ar  flame  obtained  using  the  microprobe  sampling.  The  profiles 
in  Fig.  7  are  presented  with  a  sampling  shift  Az=0.4  mm.  The  fact  that  the  value  Az  was 
determined  correctly  is  supported  by  a  satisfactory  similarity  of  temperature  (Fig.  6)  and  the 
concentration  profiles  (Fig.  7).  It  is  obvious  that  an  addition  of  TMP  increases  the  width  of  the 
zone  for  CH4,  O2,  and  H2O.  A  satisfactory  agreement  between  simulated  and  measured 
concentration  profiles  presented  in  Fig.  7  is  observed. 

The  structure  of  the  lean  CFLi/CF/Ar  flame  doped  with  2200  ppm  of  TMP  and  without 
additive  was  studied  by  MBMS.  Figure  8  shows  simulated  and  measured  concentration 
profiles  of  CH4,  O2  and  H2O  and  experimental  profile  of  temperature  in  the  lean  undoped 
flame.  Here  and  below  all  concentration  profiles  are  shifted  towards  the  burner  surface  for  the 
value  of  sampling  shift  [22].  The  sampling  shift  (z)  according  to  [22]  depends  on  the  flame 
temperature  as  Az=  zo(To/T)°‘5  where  To  -  is  initial  temperature  of  the  combustible  mixture,  T 
-  is  a  flame  temperature  for  a  certain  distance  from  the  burner,  zo  -  sampling  shift  at  T=To 
[22].  In  present  study  the  sampling  shift  varied  from  0.51-0.27  mm.  A  comparison  shows  the 
results  of  simulation  to  be  in  a  satisfactory  agreement  with  experimental  data.  This  fact,  in 
particular,  confirms  a  validity  of  this  approach  for  study  of  flame  structure  at  a  pressure  of  1 
atm. 

Figure  9  demonstrates  simulated  and  measured  concentration  profiles  of  CH4,  O2, 
H2O,  CO  and  CO2  and  measured  profile  of  temperature  in  the  lean  flame  doped  with  2200 
ppm  of  TMP.  An  agreement  between  modeling  and  experimental  data  is  observed. 

In  Figures  10  and  11  concentration  profiles  of  CH4,  O2,  and  H2O  and  measured 
profiles  of  temperature  in  the  rich  CH4/02/Ar  flame  without  additive  and  in  that  doped  with 
2200  ppm  of  TMP  are  shown  correspondingly.  As  it  is  seen  in  Figs.  10  and  11,  the  additive 
increases  the  combustion  zone  width  of  the  flame  for  about  0.4  mm.  The  increase  of 
combustion  zone  results  in  decrease  of  the  temperature  gradient  near  the  burner  surface  that  in 
its  turn  reduces  the  heat  losses  from  the  flame  into  the  burner.  That  is  why  the  inhibited 
flames  are  perturbed  much  stronger  by  a  probe  than  uninhibited  ones.  Finally  it  results  in 
appreciable  obstacles  in  reproducing  the  experimental  results  obtained  using  different  probes. 
The  use  of  a  probe  with  just  the  same  characteristics  (wall  thickness,  diameter  of  orifice)  can 
help  to  overcome  these  difficulties.  The  simulated  concentration  profiles  of  CH4,  O2,  H2O  in 
these  flames  are  presented  in  Figs.  10  and  11.  Comparison  of  experimental  and  modeling 
results  reveals  a  satisfactory  agreement. 
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Concentration  profiles  of  TMP  and  phosphorus-containing  products  of  its  destruction 
in  the  lean  CH4/O2MX  flame  doped  with  2200  ppm  PO,  PO2,  HOPO,  HOPO2  and  OP(OH)3 
are  shown  in  Fig.  12.  The  energy  of  ionizing  electrons  of  12.8  eV  for  PO,  POo  and  HOPO, 
14.5  eV  for  HOPO2  and  17.5  eV  for  OP(OH)3  made  it  possible  to  use  the  calibration 
coefficients  for  the  species  obtained  earlier  [8].  The  molecular  ion  H0P02+  and  fragmentary 
ions  of  TMP  and  some  intermediate  products  of  TMP  destruction  in  flames  such  as 
dimethylphosphate,  dimethylphosphite  and  others  contribute  in  the  peak  at  80  amu  measured 
at  energy  of  ionizing  electrons  of  14.5  eV  As  appearance  potentials  of  fragmentary  ions  of 
OPC  are  lower  than  ionization  potential  of  HOPO2  it  is  impossible  to  escape  contribution  of 
fragmentary  ions  into  peak  at  80  amu.  To  obtain  concentration  profile  of  HOPO2  it  is 
necessary  to  subtract  the  contribution  all  OPC  -  intermediates  of  TMP  destruction  in  flame  in 
intensity  of  peak  at  80  amu.  Unfortunately,  it  is  impossible  as  mass  spectra  of  some  of  the 
intermediates  (methylphosphate,  methylphosphite  and  others)  are  unknown  because  of  their 
instability  and  therefore  inaccessibility  at  normal  conditions.  Taking  into  account  the 
contributions  in  peak  at  80  amu  from  TMP  and  dimethylphosphate  (peak  at  126  amu), 
assuming  absence  of  the  contributions  into  this  peak  from  other  OPC  intermediates  and 
considering  the  concentration  profile  to  be  smooth  and  relying  on  equation  of  conservation  on 
phosphorus,  the  profiles  of  concentrations  of  HOPO2  in  the  lean  and  rich  CH4/02/Ar  flames 
were  plotted  in  Figs.  12  and  13.  So,  Figures  12  and  13  represent  final  concentration  of  HOPO2 
in  the  flames  reliably  only  at  4  mm  from  the  burner  surface  where  there  are  no  OPC. 
Concentration  profiles  of  PO,  PO2,  HOPO,  HOPO2  and  OP(OH)3  in  the  rich  CHVOo/Ar  flame 
doped  with  2200  ppm  of  TMP  are  presented  in  Fig  12.  Orthophosphoric  acid  is  an 
intermediate  product  of  TMP  destruction  in  the  rich  flame.  It  may  be  explained  by  relatively 
high  final  temperature  of  the  flame.  Unlike  the  lean  flame  in  the  rich  one  the  main  final 
phosphorus-containing  compound  is  HOPO.  One  can  notice  that  the  width  of  reaction  zone  of 
phosphorus-containing  species  is  larger  than  that  of  stable  species  (CH4,  O2,  H2O)  especially 
in  the  lean  flame.  Simulated  profiles  of  concentration  of  PO,  PO2,  HOPO,  HOPO2  and 
OP(OH)3  in  the  lean  and  rich  flames  are  presented  in  Figs.  12  and  13  correspondingly. 
Considering  difficulties  emerging  at  measuring  concentration  of  these  species  and  accuracy  of 
determination  of  their  calibration  coefficients  [8]  one  can  observe  a  satisfactory  agreement 
between  experiment  and  modeling  for  the  lean  flame.  The  structure  of  this  flame  was 
simulated  using  2  published  kinetic  models  describing  the  action  of  OPC  on  combustion.  The 
first  one  elaborated  by  Glaude  et  al.  [18,19]  is  defined  as  mechanism  2,  and  the  second  one 
proposed  by  Babushok  and  Tsang  [20,21]  is  defined  as  mechanism  3.  The  calculated  final 


16 


concentrations  of  PO,  PO2,  HOPO,  HOPO2  and  OP(OH)3  in  the  lean  and  rich  CfLi/Ch/Ar 

flames  doped  with  2200  ppm  of  TMP  are  presented  in  Tables  3  and  4. 

TABLE  3.  Fraction  of  number  of  moles  of  a  phosphorus-containing  compound  in  total 
number  of  moles  of  PCS  (in  %)  in  post- flame  zone  of  the  lean  CTL/CL/Ar  flame  calculated 
using  3  different  mechanisms,  equilibrium  state  and  measured  experimentally 


Mole  fractions 

Equilibrium 

Experimental 

Species 

1  2  3 

concentrations 

results 

at  1685  K 

Range 

PO 

0.0 

0.0 

0.4 

0.0 

0-0.6 

P02 

10.6 

0.5 

37.7 

1.0 

1-5 

HOPO 

8.3 

1.5 

7.7 

0.6 

4.8 -8.8 

HOPO2 

74.3 

96.6 

47.5 

90.6 

58  -  100 

PO(OH)3 

6.8 

0.6 

4.4 

7.8 

4.5  -  13.5 

TABLE  4.  Fraction  of  number  of  moles  of  a  phosphorus-containing  compound  in  total 
number  of  moles  of  PCS  (in  %)  in  post-flame  zone  of  the  rich  flame  CTL/CL/Ar 
(0.075/0.125/0.80)  calculated  using  3  different  mechanisms,  equilibrium  state  and 
measured  experimentally  (L=4  mm) 


Mole  fractions 

Equilibrium 

concentrations 

at  1800K 

Experimental 

results 

Range 

opecies 

1 

2bab 

3 

gloud 

PO 

4.2 

5.7 

0.0 

1.2 

1. 9-9.5 

P02 

8.1 

25.3 

6.7 

7.0 

2.2-11.0 

HOPO 

60.7 

60.3 

52.8 

82.9 

62.6-93.8 

HOPO2 

26.1 

8.3 

34.7 

8.6 

6.6-12.1 

PO(0)3 

1.0 

0.3 

0.1 

0.3 

0.0 

Final  concentrations  of  phosphorus-containing  species  in  the  flame  calculated  using 
our  model  (mechanism  1)  are  shown  in  Tables  3  and  4  together  with  experimental  data  and 
results  of  calculations  of  the  equilibrium  concentrations  of  these  species  at  final  flame 
temperature  using  thermochemic  data  accepted  in  mechanism  L.  Comparison  of  modeling 
data  with  experiment  revealed  that  all  three  models  provide  a  satisfactory  agreement  with 
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approximately  the  same  deviation  for  final  concentrations  of  phosphorus-containing  species  in 
the  lean  flame.  Comparison  of  the  results  of  thennodynamic  calculations  with  those  of 
experiment  and  modeling  demonstrated  a  distinction  between  them.  It  indicates  that 
thermodynamic  equilibrium  in  these  conditions  is  not  achieved.  Only  concentrations  of  some 
of  these  phosphorus-containing  compounds  are  close.  It  can  be  expected  considering 
overequilibrium  concentration  of  H  and  atoms  and  OH  radicals,  which  react  with  phosphorus- 
containing  species  producing  the  inhibition  effect.  Calculation  of  the  equilibrium 
compositions  involving  all  species  from  mechanism  3  gave  the  same  results.  Thermodynamic 
calculations  using  thermochemical  data  applied  in  mechanism  2  showed  that  at  the 
temperature  range  1200  K  and  higher  P(OH)3  was  the  only  phosphorus  component.  But 
modeling  and  experiment  did  not  confirm  fonnation  of  this  compound  in  post-flame  zone. 

Combustion  is  well  known  to  be  a  process  controlled  by  radicals.  Therefore,  a 
mechanism  of  action  all  inhibitors  including  OPC  on  combustion  is  based  on  effect  of 
inhibitor  and  (or)  its  destruction  products  in  flame  on  concentarion  of  active  flame  species. 
Thus,  a  study  of  influence  of  OPC  additives  on  profiles  of  temperature,  concentration  of 
stable  and  active  (H,  O,  OH)  species  in  CH4/O2  flames  and,  besides,  active  phosphorus- 
containing  destruction  products  P0,P02,  HOPO,  HOPO2  and  OP(OH)3  is  the  foreground  task 
in  understanding  the  mechanism  of  influence  of  OPC  on  combustion.  If  kinetic  model  of 
combustion  of  a  CH4/O2  doped  with  OPC  correctly  predicts  the  concentration  profiles  of 
atoms  and  free  radicals  and  phosphorus  oxides  and  acids  this  model  can  be  recognized  to  be 
perfect.  Determination  of  concentration  of  H,  O  and  OH  in  flames  is  based  on  partial 
equilibrium  of  three  "fast"  reactions.  The  reaction  and  their  equilibrium  constants  dependence 
on  temperature  [36,  37]  are  given  below: 


H2  +  OH  =H20  +  H 

K ! =0 , 1 1 3  T°’0839e7680/T 

(1) 

H2  +  0  =  H  +  OH 

K2  =  l,8T0,027e'917/T 

(2) 

02  +  H  =  OH  +  0 

K3=302T'0,374e'8620/T 

(3) 

The  detailed  equilibrium  in  flames  at  low  and  moderate  pressures  result  in  formation  of  H,  O 
and  OH  in  extra-equilibrium  concentrations  in  post-flame  zone.  The  partial  equilibrium  and 
approximation  of  quasi-stationary  concentrations  of  radicals  made  it  possible  to  express  their 
concentrations  as  functions  of  concentrations  of  H2,  O2,  H2O  and  equilibrium  constants  of  the 
reactions  (l)-(3)  in  the  flame: 


[0U]=jK2K3[02][H2];[0]= 


K,K,[H2]\Q2\ .  [OH][Q] 

[■ H20 ]  ,L  k3[o2] 
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Thus,  the  concentrations  of  H,  O  and  OH  in  the  lean  undoped  CH4/02/Ar  flame  were 
determined  and  calibration  coefficients  of  these  species  relatively  argon  were  estimated.  This 
made  it  possible  to  obtain  concentration  profiles  of  H  and  OH  in  the  lean  undoped  flame  and 
flame  doped  with  350  ppm  of  TMP  presented  in  Figs.  14  and  15.  Measurement  of  intensities 
of  peaks  at  1  and  17  amu  (H  and  OH,  correspondingly)  was  carried  out  at  ionization  energy  of 
16.5  eV.  Figures  14  and  15  demonstrate  that  introduction  of  the  inhibitor  results  in  noticeable 
decrease  of  concentrations  of  both  species.  When  the  inhibitor  was  added  the  concentration  of 
H  decreased  for  about  40  %  and  that  of  OH  decreased  in  2  times.  The  results  obtained 
demonstrate  strong  decrease  of  concentration  of  the  most  active  species  in  atmospheric 
methane  flame  when  doped  with  OPC.  It  is  noteworthy  that  earlier  [9]  it  was  reported  about 
increase  of  OH  concentration  for  20  %  in  Hj/Oi/Ar  flame  stabilized  at  low  pressure  when  an 
additive  of  TMP  was  introduced.  But  this  phenomenon  was  observed  for  low  pressure 
hydrogen  flames  promoted  (not  inhibited)  by  additives  of  TMP.  As  to  oxygen  atom,  it  is 
impossible  to  measure  in  methane  flame  the  intensity  of  peak  at  16  amu  corresponding  only  O 
as  its  ionization  potential  (13.618  eV)  is  higher  than  ionization  potential  of  CH4  (12.98  eV) 
[38].  In  Figure  16  the  intensity  profiles  of  peak  at  15  amu  corresponding  to  CH3  in  the  lean 
C^/CF/Ar  flame  without  additive  and  in  the  same  flame  but  doped  with  350  ppm  of  TMP  are 
presented.  The  measurement  of  the  mass  peak  intensity  was  perfonned  at  ionization  energy  of 
14.2  eV.  The  profiles  presented  in  Fig.  16  show  that  the  additive  of  TMP  effects  the 
concentration  of  CH3  similarly  to  those  of  H  and  OH  decreasing  it  twice.  So,  the  data  shown 
above  indicate  that  effectiveness  of  OPC  as  inhibitors  consists  in  strong  influence  of  OPC 
concentration  of  the  most  active  flame  species. 

Structure  of  C3H.fi/O2/Ar  Flames 

Figure  17  shows  the  simulated  and  measured  profiles  of  concentration  of  stable  species  C3H8, 
O2,  H2O  and  CO2  and  the  profile  of  temperature  in  the  lean  CFFF/CF/Ar  flame  without 
additive.  Figure  18  presents  the  temperature  and  concentration  profiles  obtained  in  experiment 
and  by  modeling  for  the  lean  CFFFAVAr  flame  doped  with  1200  ppm  of  TMP.  Figures  19 
and  20  present  the  corresponding  data  for  the  rich  CFFFAF/Ar  flame  without  additive  and 
doped  with  1200  ppm  of  TMP.  Inhibition  effect  of  the  additive  results  in  appreciable  moving 
of  the  flame  front  away  from  the  burner  accompanied  with  decreases  of  heat  losses  in  to  the 
burner  as  it  was  demonstrated  for  C^/CF/Ar  flames  above.  Maximal  temperature  of  undoped 
flames  ranges  1500-1600K  whereas  the  temperature  of  the  doped  flames  ranges  1700-1800  K. 
The  temperature  rise  is  explained  by  decrease  of  heat  losses  into  the  burner  in  the  doped 
flames  due  to  increasing  of  combustion  zone  rather  than  due  to  effect  of  OPC  on 
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thermophysical  properties  of  the  flame.  A  kinetic  model  for  destruction  of  TMP  together  with 
Konnov  mechanism  of  propane  oxidation  provides  an  appropriate  agreement  between 
measured  and  simulated  profiles  of  concentration  of  stable  species.  Besides  concentration 
profiles  of  H  and  OH  have  been  detennined  in  lean  undoped  flame  (Fig.21).  Results  of 
modeling  calculation  of  concentration  profiles  of  H  and  OH  in  the  lean  flame  also  shown  in 
Fig.  23  demonstrate  an  agreement  with  experimental  data. 

The  measured  and  calculated  concentration  profiles  of  final  phosphorus-containing 
products  of  TMP  destruction  in  the  lean  flame  -  PO,  PO2,  HOPO,  HOPO2  and  OP(OH)3  are 
shown  in  Fig.  22.  Similar  results  for  the  rich  C3H8.O2MX  flame  are  plotted  in  Fig.  23. 
Comparing  the  data  obtained  for  both  C3H8/O2MX  flames  doped  with  TMP  one  can  observe 
appreciable  differences  between  composition  of  phosphorus-containing  species  in  the  lean 
and  rich  flames.  In  the  lean  CFFRTR/Ar  flames  the  major  final  product  of  TMP  destruction  is 
HOPO2  whereas  in  the  rich  flame  apart  from  HOPO2  such  compounds  as  HOPO,  PO2  and 
OP(OH)3  present  in  higher  concentrations.  The  difference  in  the  final  temperature  in  all 
flames  is  negligible.  Thus,  a  difference  in  concentration  of  active  flame  components,  resulting 
from  different  composition  of  combustible  mixtures,  effects  on  composition  of  phosphorus- 
containing  compounds  in  flames. 

Comparing  calculated  and  measured  concentration  profiles  of  phosphorus-containing 
species  one  can  observe  a  satisfactory  agreement  between  experiment  and  modeling  for  the 
lean  CFFR/Oi/Ar  flame  doped  with  1200  ppm  of  TMP.  As  to  the  rich  flame  an  appreciable 
disagreement  between  simulated  and  measured  concentration  profiles  is  observed.  Final 
concentrations  of  phosphorus-containing  compounds  in  the  lean  and  rich  ChFR/CA/Ar  flames 
obtained  by  experimentally  and  by  modeling  are  compared  in  Tables  5  and  6  correspondingly. 

TABLE  5.  Fraction  of  number  of  moles  of  a  phosphorus-containing  compound  in  total 
number  of  moles  of  PCS  (in  %)  in  post-flame  zone  of  the  lean  flame  calculated  using 
3  different  mechanisms,  equilibrium  state  and  measured  experimentally 


Mole  fractions  Equilibrium  Experimental 


Species 

1 

2 

3 

concentrations 

at  1760K 

results 

Range 

PO 

0.0 

0.0 

0.0 

0.0 

0-0.2 

P02 

15.5 

1.2 

29.6 

3.0 

3  -  15 

HOPO 

8.3 

2.5 

8.6 

1.9 

3.6 -6.8 

HOPO2 

73.0 

96.0 

59.4 

91.3 

60-100 

20 


PO(OH)3  3.2  0.3  2.4  3.8 


The  simulated  concentrations  tabulated  above  are  obtained  using  3  different  models:  our 
model  (mechanism  1),  the  model  of  Glaude  et  al.  [18,19]  (mechanism  2)  and  the  model  of 
Babushok-Tsang  [20,21]  (mechanism  3).  Besides,  both  tables  present  the  equilibrium 
concentration  of  the  species  at  final  flame  temperature.  It  should  be  specially  pointed  out  that 
equilibrium  concentrations  were  calculated  using  thermochemical  data  recommended  by 
Melius  [4],  which  is  used  in  our  kinetic  mechanism  and  the  mechanism  of  Babushok-Tsang 
[20,21]. 

TABLE  6.  Fraction  of  number  of  moles  of  a  phosphorus-containing  compound  in  total 
number  of  moles  of  PCS  (in  %)  in  post-flame  zone  of  the  rich  flame  calculated  using  3 
different  mechanisms,  equilibrium  state  and  measured  experimentally 


Species 

Mole  fractions 

1  2  3 

Equilibrium 

concentrations 

at  1700  K 

Experimental 

results 

Range 

PO 

4.3 

0.0 

0.2 

0.8 

0-1.1 

P02 

6.6 

4.9 

14.1 

4.8 

2  -  11.5 

HOPO 

63.9 

61.6 

75.5 

85.2 

22-41 

H0P02 

23.8 

32.7 

9.6 

8.7 

34-63 

PO(OH)3 

1.4 

0.1 

0.6 

0.5 

5  -  15 

The  application  of  thermochemistry  developed  by  Glaude  and  others  predicts  all  phosphorus 
to  transfonn  into  P(OH)3  that  is  not  confirmed  by  modeling.  The  data  presented  in  Tables  2 
and  3  show  that  all  the  mechanisms  provide  approximately  the  same  agreement  of  modeling 
data  with  experimental  results  for  the  lean  flame,  whereas  none  of  the  mechanisms  is  able  to 
predict  the  concentrations  of  phosphorus-containing  species  with  an  appropriate  accuracy. 
Tables  2,  3  also  revealed  that  equilibrium  concentrations  of  these  species  are  not  achieved  in 
both  flames.  Only  concentration  of  some  species  is  close  to  equilibrium  ones. 

Measurements  of  inhibition  effect  of  OPC  on  burning  velocity 

Study  the  effect  of  number  of  OPC  including  a  number  of  novel  substances  on  burning 
velocity  of  stoichiometric  CPfi/air  and  C3Hx/air  flames  have  been  carried  out  at  atmospheric 
pressure.  TMP  loading  dependencies  of  burning  velocity  of  a  stoichiometric  CPfi/air  flame 
measured  experimentally  and  calculated  using  our  kinetic  model  and  4  different  mechanisms 
of  methane  oxidation  (see  Table  2)  are  presented  in  Fig.  24.  Modeling  calculations  showed 
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CHEMKIN  mechanism  to  provide  a  better  agreement  with  experimental  results.  Figure  25 
demonstrates  TMP  loading  dependencies  of  burning  velocity  of  a  stoichiometric  CH4/air 
flame  measured  experimentally  and  calculated  using  GRI  3.0  mechanism  and  3  published 
kinetic  models  for  action  of  phosphorus-containing  species  on  combustion.  As  one  can  see  all 
3  models  involving  phosphorus  substances  give  very  close  results.  Analyzing  Figs.  24  and  25 
one  can  conclude  that  a  disagreement  of  modeling  results  with  experimental  data  can  be 
explained  possibly  by  drawbacks  of  GRI  3.0  rather  than  the  phosphorus-involving  models. 
Nevertheless,  a  disagreement  with  experiment  can  be  decreased  by  increasing  the  pre¬ 
exponential  factor  of  reaction  H+P02+M=H0P0+M  in  4  times.  The  modeling  results 
obtained  using  mechanism  1  with  modified  rate  constant  (defined  as  mechanism  1.1)  is 
presented  in  Fig.  25. 

TMP  loading  dependence  of  nonnalized  burning  velocity  of  a  stoichiometric  C3fU/air 
flame  is  shown  in  Fig.  26.  Modeling  using  3  published  mechanisms  involving  OPC 
destruction  (Fig.  26)  is  in  a  satisfactory  agreement  with  experiment.  Comparing  TMP  loading 
dependence  of  burning  velocity  of  CfU/air  and  CUHg/air  stoichiometric  flames  (Fig.  27)  one 
can  see  that  TMP  loading  dependence  of  the  burning  velocity  does  not  differ  drastically  for 
both  methane  and  propane  flames.  Nevertheless,  an  inhibition  effect  of  TMP  on  CfU/air  flame 
is  slightly  stronger. 

The  results  of  measurement  of  burning  velocity  of  CfU/air  flames  of  various 
equivalence  ratio  (in  the  range  of  measurement  of  cp  (0.9- 1.5)  for  the  flame  without  additive 
and  doped  with  520  ppm  of  TMP  (Fig.  28).  Burning  velocities  were  normalized  on  the 
maximum  value  of  burning  velocity  of  the  undoped  flame.  In  Fig.  28  the  value  f=(uundoped- 
UdoPed)/uUndoped,  which  presents  a  relative  decrease  of  burning  velocity  when  520  ppm  of  TMP 
is  added  to  the  flame  and  characterizes  effectiveness  of  inhibitor  action  on  a  flame. 
Equivalence  ratio  dependence  of  f  passes  though  maximum  at  (p=  1.1 -1.2 

The  fluorine-  and  nitrogen-containing  volatile  (in  comparison  with  regular  phosphates 
and  phosphonates)  compounds  are  of  advanced  interest.  The  studied  OPC  and  some  of  their 
properties  are  presented  in  Table  7. 


Table  7.  Organophosphorus  compounds,  which  were  tested  as  flame  inhibitors. 


Compound 

Boiling  point  [UC]  at  pressure  [Torr] 

(C2H50)3P0 

215/760 

(CF3CH20)3P0 

73/8 

[(CH3)2N]3PO 

123/19 
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(CH30)3P0 


180/760 


Dependencies  of  inhibiting  effectiveness  of  OPC  tabulated  above  for  propane  flame  are 
presented  in  Fig.  29.  The  data  shown  in  Fig.  29  demonstrate  that  all  tested  OPC  have 
approximately  the  same  flame  inhibiting  effectiveness.  F  and  N  atoms  introduced  in 
molecules  of  some  OPC  do  not  influence  on  properties  of  OPC  as  flame  inhibitors. 


Measurement  of  extinction  strain  rate  in  a  non-premixed  CH_/02  flame 


The  extinction  strain  rate  in  a  non-premixed  methane  flame  with  OPC  loading  stabilized  on 
the  opposed-jet  burner  at  1  atmosphere  was  measured.  The  opposed-jet  burner  was  designed 
in  accordance  with  recommendations  of  [39].  The  internal  diameter  of  nozzles  was  6.8  mm, 
the  distance  between  nozzles  was  6.8  mm.  The  flame  stabilized  on  the  fuel  side  of  stagnation 
plane.  The  sheath  flow  was  Ar.  The  maximum  of  flow  rate  was  50  cm  /s  for  both  fuel  and 
oxidizer.  The  burner  was  aligned  upright.  Flow  rates  of  the  components  are  presented  in  Table 
8. 


Table  .8.  Volumetric  flow  rate  (NTPcm 7s)  of  the  components  of  combustible  mixture. 


FUEL 

OXIDIZER 

ch4 

n2 

o2 

n2 

10 

40 

20 

30 

The  strain  rate  was  calculated  according  to  equation  [40]: 

£  {'  K,  1  pm  J 

here  V  is  the  stream  velocity  and  p  is  the  stream  density;  L  is  the  separation  distance  between 
nozzles.  The  typical  value  of  strain  rate  for  flames  varies  from  100  to  800  s"1  [39,41,42]  and 
depend  from  the  composition  of  fuel  and  oxidizer  and  burner  design.  In  our  experiment  the 
value  of  strain  rate  without  loading  OPC  was  750  s'1.  The  extinction  strain  rates  were 
measured  by  increasing  the  flow  rates  of  fuel  and  oxidizer  until  flame  is  quenched.  The 
loading  of  TMP  ((CH30)3P0)  was  from  500  to  3650  ppm  and  TEP  (((C2H50)3P0)  was  from 
500  to  1200  ppm.  In  Fig.  30  nonnalized  extinction  strain  rate  of  the  flame  without  additive 
and  doped  with  of  OPC  are  shown  as  a  function  of  OPC  loading.  The  data  obtained 
demonstrate  that  TMP  and  TEP  loading  dependence  of  the  burning  velocity  practically  are 
equal  for  both  compounds. 


2F 


1  -l- 


'fuel 
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CONCLUSION 

Sensitivity  analysis  of  phosphorus-involving  reactions  of  our  mechanism  was 

performed  to  determine  the  reactions,  which  control  the  burning  velocity  of  a  doped  CHVair 

flame.  Six  reactions  are  specified  in  Fig.  31.  The  stage  H+P02+M=H0P0+M  was 

demonstrated  to  be  the  most  important  (sensitivity  coefficient  -35).  The  rate  constants  of  this 

reaction  accepted  by  the  authors  of  different  kinetic  mechanisms  are  tabulated  below. 

Table  9.  Rate  constant  (expressed  as  ATne'E/kT)  of  the  reaction  H+P02+M=H0P0+M  in 
different  kinetic  models 


MECHANISM 

A 

n 

E  (kal/mol) 

A/ Ajwarowski 

our  mechanism 

9.73X1024 

-2.0 

645.0 

10.0 

Glaude  et  al.  [18,19] 

1.46xl025 

-2.0 

645.0 

15.0 

Babushok-Tsang  [20,21] 

5.00X1024 

-2.0 

645.0 

5.1 

Twarowski  [5-7] 

9.73X1023 

-2.0 

645.0 

1.0 

The  data  presented  in  Table  9  indicate  that  pre-exponential  factor  of  rate  constant  of  this 
reaction  accepted  in  all  models  except  Twarowski  one  are  very  close  (differ  within  2.25 
times).  Besides  the  temperature  dependence  of  the  rate  constants  is  the  same  in  all  the  models. 
That  is  why  the  results  of  calculation  of  burning  velocity  using  these  mechanisms  are  very 
close.  It  means  that  a  essentially  different  models  but  having  a  close  rate  constant  of  the  key 
reactions  can  provide  a  good  agreement  with  experiment.  It  is  noteworthy  that  all  authors  do 
not  accept  the  value  of  rate  constant  of  this  reaction  recommended  by  Twarowski.  Another 
important  reactions  and  their  rate  constants  recommended  by  different  authors  are  tabulated 
below. 

Table  10.  The  most  important  reactions  responsible  for  inhibition  effect  and  their  rate 
constants  in  our  model,  model  of  Glaude  et  ah,  model  of  Babushok-Tsang 


expressed  as  k  =  A  Tn*exp(-E/RT)) 


Reaction 

our  model 

Glau 

Babus 

de  et  al. 

hok-Tsang 

A 

E 

A 

A 

0H+P02+M=H0P02+M 

1.6xl025 

-2.28 

285.0 

1.6xl024 

1.6xl024 

h+hopo2=h2o+po2 

6.32X1011 

0.00 

11930.0 

6.32X1012 

6.32X1012 

h+po2+m=hopo+m 

9.73X1024 

-2.04 

645.0 

9.73X1025 

9.73X1024 

oh+hopo=h2o+po2 

3.16X1011 

0.00 

0.0 

3.16X1012 

3.16X1012 

h+hopo=h2+po2 

7.9xlOn 

0.00 

43.0 

7.9xl012 

7.9xlOn 

o+hopo=oh+po2 

1.58X1012 

0.00 

0.0 

1.58xl013 

1.58X1013 

o+hopo+m=hopo2+m 

1.3xl024 

-2.11 

995.0 

1.3xl023 

1.3xl023 

o+hopo2=o2+hopo 

6.32X1011 

0.00 

8236.0 

6.32X1012 

6.32X1012 

Units  are  mol,  cm3  ,  s,  K  and  calxmol 
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Temperature  dependence  for  all  rate  constants  is  recommended  by  Twarowski  and  the 
same  in  all  3  models.  To  analyze  if  the  key  stages  are  the  same  for  rich  and  lean  C^ffi/air 
flames  sensitivity  coefficients  on  burning  velocity  were  calculated  (Table  1 1). 


Table  11.  The  coefficients  of  sensitivity  for  burning  velocity 
(in  %)  of  the  lean  (cp=0.8)  and  rich  (cp=  1 .2)  CstR/air  flames. 


N 

Reaction 

c,=(u0-uk/5/u)  xl00% 

(p=0.8 

(p=1.2 

1. 

h+po2+m=hopo+m 

-12.7 

-12.2 

2. 

o+hopo=po2+oh 

-4.7 

-4.2 

3. 

oh+po2+m=hopo2+m 

3.0 

0.8 

4. 

oh+hopo=h2o+po2 

-1.6 

-4.9 

5. 

h+hopo2=h2o+po2 

-1.2 

-3.2 

6. 

o+hopo2=o2+hopo 

-1.1 

2.0 

Where  uo  -  burning  velocity  at  recommended  rate  constant;  Uk/5  -  burning  velocity  calculated 
with  the  rate  constant  reduced  in  5  times.  Table  1 1  shows  that  the  sensitivity  coefficients  of  a 
reaction  on  burning  velocity  may  depend  on  equivalence  ratio  of  a  flame.  Thus,  reaction  4 
(Table  1 1)  is  the  second  one  on  importance  for  the  rich  CoFR/air  flame  and  the  fourth  for  the 
lean  flame. 

As  phosphorus-containing  products  of  OPC  destruction  in  a  flame  are  known  to  be 
responsible  for  inhibition  effect  an  appreciable  efforts  were  made  to  measure  and  simulate 
their  concentration  profiles.  Two  other  published  models  for  OPC  destruction  in  a  flame  were 
validated  by  comparing  simulated  and  experimental  profile  of  concentration  of  these  products. 
There  are  several  conclusions,  which  can  be  done. 

None  of  the  tested  models  can  predict  concentration  profiles  of  PO,  PO2,  HOPO,  HOPO2  and 
OP(OH)3  in  the  rich  flames.  It  can  be  connected  with  insufficient  understanding  of 
combustion  chemistry  of  rich  flames  as  well  as  combustion  chemistry  of  OPC.  As  an 
illustration  several  examples  can  be  cited.  Possibility  of  presence  of  the  species,  which  are  the 
products  of  addition  of  hydrocarbon  radicals  to  phosphorus  oxides  in  rich  flames  was  reported 
in  paper  [34].  A  variety  of  PCS  of  complex  structure  was  found  in  combustion  products  of 
rich  C2H2/02/Ar  flames  doped  with  TMP  and  tri(dimethylamino)phosphine.  Results  reported 
in  [35]  indicate  that  combustion  chemistry  of  OPC  is  not  well  understood.  Thus,  the 
experimental  data  obtained  demonstrated  that  an  additive  of  DMMP  to  the  lean  and  rich 
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methane-oxygen  atmospheric  flame  results  in  rise  of  concentration  of  CH3OH  in  10  times. 
But  results  of  modeling  using  mechanisms  2  and  3  failed  to  describe  this  fact. 

The  models  for  oxidation  of  hydrocarbons,  especially  those  involving  C2-C4  species 
may  contain  rate  constants,  which  are  determined  insufficiently  accurate.  The  influence  of  a 
model  for  methane  oxidation  on  final  concentration  of  phosphorus-containing  species  in  the 
lean  and  rich  CTLi/Cb/Ar  flame  doped  with  2200  ppm  of  TMP  can  be  demonstrated  by  the 
following  data  shown  in  Tables  12  and  13. 

Table  12.  Final  concentrations  of  phosphorus-containing  species  in  lean  CFB/Ch/Ar  flame 
doped  with  2200  ppm  of  TMP  calculated  using  our  inhibition  mechanism  and  3  different 
mechanisms  of  CH4  oxidation. 


Species 

Experim 

ental  data 

Simulated  final  concentrations  using  our 
model  for  TMP  destruction  and  3  different 
models  for  CH4  oxidation 

CHEMKIN 

Konnov 

GRI  3.0 

PO 

0-0.6 

0 

0 

0 

P02 

1-5 

11.7 

8.7 

7.7 

HOPO 

4.8 -8.8 

4.4 

3.1 

2.7 

HOPO2 

58-100 

76.8 

81.5 

82.0 

H3PO4 

4.5  -  13.5 

7.1 

7.6 

7.6 

Table  13.  Final  concentrations  of  phosphorus-containing  species  in  rich  CFB/Ch/Ar  flame 
doped  with  2200  ppm  of  TMP  calculated  using  our  inhibition  mechanism  and  3  different 
mechanisms  of  CH4  oxidation. 


species 

Experim 

ental  data 

Final  concentrations  obtained  using  our  model 
for  TMP  destruction  and  3  models  for  CH4 
oxidation 

CHEMKIN 

Konnov 

GRI  3.0 

PO 

1. 9-9.5 

4.2 

2.6 

2.9 

P02 

2.2-11.0 

8.1 

7.9 

7.8 

HOPO 

62.6-93.8 

60.7 

50.7 

53.3 

HOPO2 

6.6-12.1 

26.1 

37.3 

34.6 

H3PO4 

0 

1.0 

1.5 

1.3 

Although  the  concentrations  simulated  using  different  methane  oxidation  mechanisms  differ 

not  more  than  in  1.5-2  this  fact  is  of  interest.  The  modeling  results  can  be  explained  by 

different  concentrations  of  active  species  predicted  by  different  mechanisms  of  methane 

combustion  as  phosphorus-containing  components  of  a  flame  react  actively  with  free  atoms 

and  radicals.  As  the  concentration  of  active  species  in  the  studied  flames  is  higher  than 

summary  concentrations  of  phosphorus  species  a  small  change  in  concentration  can  result  in 

appreciable  change  of  concentrations  of  phosphorus-containing  products. 

Table  14.  Calculated  maximal  concentration  of  active  species  in  the  in  the  rich  CFB/Ch/Ar 
flame  doped  with  2200  ppm  TMP. 

"Species  CHEMKIN  Konnov  GRI3.0 

H  1.61xl0'3  9.75xl0'4  1.03xl0'3 
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OH  4.64x1  O'4  3. 94x1  O'4  3.87xl0'4 

_0 _ 8.34xl0'5 _ 4.32xl0'5  5.45xl0'5 

Thus,  it  is  difficult  to  expect  an  agreement  between  simulated  and  measured  concentration  of 

phosphorus-containing  species  better  than  ±50  %  for  several  flames  simultaneously.  In  spite 

of  necessity  of  further  modification  of  the  models  describing  inhibition  action  of  OPC  on  a 

flame,  the  problems  with  a  kinetic  mechanism  of  hydrocarbon  oxidation  should  be  taken  into 

consideration. 

The  ways  of  further  work  on  refining  of  the  models  for  inhibition  of  combustion 
should  be  discussed.  The  updating  a  model  with  new  stages  involving  phosphorus-containing 
compounds  seems  to  be  a  good  idea.  But  the  example  of  published  models  demonstrated  that 
an  introduction  of  new  stages  describing  the  interaction  of  phosphorus  oxides  with  carbon- 
containing  species  was  not  quite  helpful.  The  introduction  of  additional  stages  should  be  well- 
founded  and  directed  to  solution  of  definite  problems  and  based  on  sensitivity  analysis. 
Another  way  of  refinement  consists  in  calculation  of  rate  constants  of  reactions.  Although  this 
way  is  very  difficult  and  consumes  a  lot  of  time  it  seems  to  be  on  of  prospective. 
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SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

♦  Methodic  difficulties  of  MBMS  measurement  of  the  structure  of  atmospheric 
flames  were  analyzed  and  recommendations  were  elaborated.  The  method  of 
evaluation  of  perturbation  of  an  atmospheric  flame  caused  by  a  sonic  probe  was 
proposed. 

♦  Profiles  of  temperature  and  concentration  of  following  species  in  the  lean  (cp=0.8) 
and  rich  (cp=1.2)  CH4/02/Ar  and  C3H8/02/Ar  flames  were  measured:  CH4,  C3H8, 
02,  H20,  CO,  C02,  H,  OH,  TMP,  PO,  P02,  HOPO,  H0P02  and  OP(OH)3. 

♦  TMP  loading  dependence  of  burning  velocity  of  stoichiometric  CH4/air  and 
C3H8/air  flames  were  measured  using  Mache-Hebra  nozzle  burner  and  the  total 
area  method  of  a  flame  image. 

♦  TMP  and  triethylphosphate  loading  dependence  of  normalized  extinction  strain  rate 
of  opposed-jet  CH4/02/N2  flame  was  measured.  Both  OPC  demonstrated  close  fire 
suppression  effectiveness. 

♦  The  structure  of  the  lean  ((p=0.8)  and  rich  ((p=1.2)  CH4/02/Ar  and  C3H8/02/Ar 
flames  without  additives  and  doped  with  TMP  and  burning  velocity  of  methane-air 
and  propane-air  flames  at  various  concentration  of  TMP  additives  were  simulated 
using  PREMIX  and  CHEMKIN  codes. 

♦  Earlier  proposed  kinetic  mechanism  for  TMP  destruction  in  low  pressure  H2/02/Ar 
and  CH4/02/Ar  flames  was  refined  basing  on  comparison  of  measured  and 
calculated  results  on  the  structure  and  burning  velocity  of  atmospheric  flames. 

♦  Two  published  models  for  destruction  of  OPC  in  a  flame  proposed  by  Glaude  et  al. 
and  Babushok  and  Tsang  were  used  for  simulation  of  the  structure  and  burning 
velocity  of  the  flames. 

♦  The  sensitivity  analysis  was  applied  to  specify  the  stages  responsible  for  inhibition 
and  their  rate  constants  accepted  for  all  3  models  have  been  compared. 

♦  The  mechanisms  for  methane  combustion  available  in  literature  were  used  for 
simulation  of  concentration  profiles  of  phosphorus-containing  species  and  burning 
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velocity  of  a  stoichiometric  CH4/air  flame.  Their  influence  on  composition  of 
phosphorus-containing  flame  products  was  demonstrated. 

♦  Recommendations  for  further  investigation  of  combustion  chemistry  of  OPC  and 
refinement  of  the  kinetic  model  for  OPC  destruction  in  flame  have  been  proposed. 
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APPENDIX 

Figures 


Fig.l.  Temperature  profiles  in  lean  CF^/CVAr  (0.06/0.15/0.79)  flame  at  1 
bar.  Curve  1 -without  probe(S=°o),  curve  2-8=0. 3mm,  curve  3-8=0. 13mm. 


Fig.2.  Temperature  profiles  in  lean  flame  CH4/02/Ar  (0.06/0.15/0.79)  doped 
with  2200ppm  of  TMP  at  1  bar.  Curve  1 -without  probe  (S=°°),  curve  2- 
8=0. 25mm,  curve  3-Al=0.1mm. 


Fig. 3.  Temperature  profiles  in  rich  flame  CH4/02/Ar  (0.075/0.125/0.8)  at  1  bar. 
Curve  1 -without  probe(8=oo),  curve  2-8=0. 3mm. 
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Fig.  4.  Temperature  profiles  in  the  lean  CH4/02/Ar  flame  doped  with  1200 
ppm  of  TMP  measured  by  a  thermocouple  placed  at  distances  8=0.6  and 
0.7  mm  from  the  probes  and  by  pneumatic  probes. 
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Fig.  5.  Normalized  profiles  of  temperature  aT=(T-T0)/(TrTo)  and 
concentration  of  02  0CC= ( C  -  C  0 )/( C  r  C  0 )  in  a  near  stoichiometric  CH4/02/Ar 
flame  without  additive  of  TMP.  Line  -  temperature;  filled  symbols  - 
results  of  modeling;  open  symbols  -  results  of  experiment. 


Fig.  6.  Temperature  profiles  in  the  lean  unperturbed  CH4/02/Ar  flame  without 
additive  (open  symbols)  and  doped  with  350  ppm  of  TMP  (filled  symbols). 


Fig.  7.  Concentration  profiles  of  CH4,  O2,  and  FFO  in  the  lean 
CH4/02/Ar  flame  without  additive  (open  symbols)  and  doped  with  350 
ppm  of  TMP  (filled  symbols)  obtained  by  microprobe  technique;  lines  - 
modeling. 


Fig. 8  The  concentration  profiles  for  undoped  lean  CH4/02/Ar 
(0.06/0.15/0.79)  flame.  Solid  lines  -  modeling;  symbols  -  experiment 
(MBMS  technique). 
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Fig. 9  The  concentration  profiles  for  doped  lean  CH4/02/Ar 
(0.06/0.15/0.79)  flame  with  2200  ppm  of  TMP.  Solid  lines  - 
modeling;  symbols  -  experiment. 


Fig.  10  The  concentration  profiles  for  undoped  rich  CH4/02/Ar 
(0.075/0.125/0.80)  flame.  Solid  lines  -  modeling;  symbols  - 
experiment. 
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Fig.  11.  The  concentration  profiles  for  doped  rich  CHVCb/Ar 
(0.075/0.125/0.80)  flame  doped  with  2200  ppm  of  TMP.  Solid  lines 
-  modeling;  symbols  -  experiment  (MBMS  technique). 


Fig.  12.  Profiles  of  Concentration  of  TMP  and  phosphorus-containing  products  of 
its  destruction  in  the  lean  CH4/02/Ar  flame  doped  with  2200  ppm  of  TMP; 
symbols  -  experiment,  lines-  modeling. 
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Fig.  13  Concentration  profiles  of  TMP  and  phosphorus  containing 
products  of  its  destruction  in  the  rich  CH4/02/Ar  (0.075/0.125/0.80) 
and  sum  of  experimentally  measured  concentration  of  TMP, 
H3PO4,  HOP02,  HOPO,  P02,  PO.  Solid  lines  -  modeling;  symbols 
-  experiment. 
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Fig.  14.  Concentration  profiles  of  H  and  OH  in  the  lean  CH4/02/Ar 
flame  without  additive. 
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Fig.  15.  Concentration  profiles  of  H  and  OH  in  the  lean 
CH4/02/Ar  (0.06/0.15/0.79)  flame  doped  with  350  ppm  of  TMP. 
Solid  lines  -  modeling;  symbols  -  experiment. 
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Fig.  16.  Profiles  of  intensities  of  peak  at  15  amu  (CH3)  in  the  lean 
CH4/02/Ar  (0.06/0.15/0.79)  flame  without  additive  (open  symbols) 
and  doped  with  350  ppm  of  TMP  (filled  symbols). 
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Fig.  17  Profiles  of  temperature  (experiment)  and  concentration  of  stable  species 
the  lean  QHyCb/Ar  flame  without  additive.  Symbols  -  experiment,  lines 
modeling. 
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Fig.  18  Profiles  of  temperature  (experiment)  and  concentration  of  stable  species  in 
the  lean  CsFlg/CF/Ar  flame  doped  with  1200  ppm  of  TMP.  Symbols  -  experiment, 
lines  -  modeling. 
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Fig.  19  Profiles  of  temperature  (experiment)  and  concentration  of  stable  species 
in  the  rich  C3H8/02/Ar  flame  without  additive.  Symbols  -  experiment,  lines  - 
modeling. 
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Fig.  20.  Profiles  of  temperature  (experiment)  and  concentration  of  stable  species 
in  the  rich  C^FIg/Cb/Ar  flame  doped  with  1200  ppm  of  TMP.  Symbols  - 
experiment,  lines  -  modeling. 


Fig.  21.  Concentration  profiles  of  ft  atoms  and  OH  radicals  in  the  lean 
C3H8/02/Ar  flame  without  additive.  Symbols  -  experiment,  lines  - 
modeling. 
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Fig.  22.  Concentration  profiles  of  TMP  and  of  final  phosphorus-containing 
products  of  its  destruction  PO,  P02,  HOPO,  FIOPCF  in  the  lean  C3H8/02/Ar  flame 


Fig.  23.  Concentration  profiles  of  final  phosphorus-containing  products  of 
TMP  destruction  and  TMP  in  the  rich  C3H8/02/Ar  (2.92/12.08/85)  flame  doped 
with  1200  ppm  of  TMP.  P=1  bar.  Solid  lines  -  modeling;  symbols  - 
experiment. 


Fig.  24.  TMP  loading  dependence  of  burning  velocity  of  CH4/air  flame 
measured  experimentally  (symbols)  and  calculated  by  modeling  (lines) 
using  our  TMP  destruction  model  and  4  various  models  of  CH4 
oxidation. 
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Fig.  25.  TMP  loading  dependence  of  burning  velocity  of  a  stoichiometric 
CH4/air.  Symbols  -  experiment;  lines  -  modeling  performed  using  GRI  3.0 
mechanism  and  mechanisms  1 ,2,  3  and  mechanism  1 . 1  with  modified  rate 
constant  H+P02+M=H0P0+M,  A=4.0xl025. 


Fig.  26.  Burning  velocity  of  the  stoichiometric  C3F18  +  air  mixture  as  a 
function  of  TMP  loading.  Symbols  -  experiment;  lines  -  modeling  results 
obtained  using  three  different  mechanisms  for  C3H8  +  air  mixture. 
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Fig.  27.  TMP  loading  dependence  of  burning  velocity  of  CH4/air 
(circle)  and  C3H8/air  (square)  flames  (experiment). 
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Fig.  28.  Equivalence  ratio  dependence  of  burning  velocity  of  CH4/air 
flames  without  additive  and  doped  with  520  ppm  of  TMP. 
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Fig.  29  Normalized  burning  velocity  of  stoichiometric  C3F18  -  Air  flame  as 
function  inhibitor  loading. 
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Fig. 30  Normalized  extinction  strain  rate  of  the  flame  without  additive 
and  doped  with  of  OPC  as  a  function  of  OPC  loading. 
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Fig.  31.  Influence  of  varied  rate  constants  of  elementary  reactions  on  calculated  value 
q=  1 00%x(u-u(k/5))/u  for  our  mechanism,  where  u-  burning  velocity  in  stoichiometric  CH 4/air 
mixture  doped  with  TMP  of  different  concentrations  at  1  bar,  u(k/5)  -  burning  velocity  for  the 
rate  constants  decreased  in  5  times. 


